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Many extensions of the standard model (SM) gauge 
group predict the existence of electrically-neutral, mas- 
sive gauge bosons commonly referred to as Z' 0, Q, y, S 
5]. The leptonic decays Z' — > £ + £~ provide the most dis- 
tinctive signature for observing the Z' signal at a hadron 
collider. In two recent publications, the Collider De- 
tector at Fermilab (CDF) Collaboration has set limits 



on different Z' models by analyzing the invariant mass 
(Mu) spectrum of the dielectron, dimuon, and ditau fi- 
nal states, using a dataset corresponding to an integrated 
luminosity of approximately 0.2 fb _1 0,0- Besides the 
dilepton mass Mu, it has been shown that the angu- 
lar distribution of the dilepton events can also be used 
to test the presence of a Z' boson by detecting its in- 
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terference with the SM Z boson In this Letter, for 
the first time at a hadron collider, the massive resonance 
search technique (Mu analysis) is extended to include 
dilepton angular information to detect Z' — > e + e~ de- 
cays in 0.45 fb _1 of data accumulated with the CDF II 
detector. As the integrated luminosity increases, the sen- 
sitivity of the standard Mu analysis plateaus; adding the 
angular information starts to become an important han- 
dle for extending the Z' exclusion reach and discovery 
potential. In this Letter many of the theoretical Z' mod- 
els are surveyed, and results are reported for the sequen- 
tial Z' SM , the canonical superstring- inspired E$ models 
Z x , Z^jZjv Zi, Z N , Z sec [lEii, the "littlest" Higgs Z H 
model [4J, |5j , the four generic model- lines of Ref. |2| , and 
contact-interaction searches. No significant evidence of 
a Z' signal is found, and the tightest constraints to date 
are set on most of these models. 

The CDF detector is described in detail elsewhere ^lj . 
For this study, the relevant subdetector systems are the 
central tracking chamber (COT) and the central and the 
plug calorimeters. The COT is a 96-layer open-cell drift 
chamber immersed in a 1.4 Tesla magnetic field, used to 
measure charged particle momenta within the pseudo- 
rapidity range \r\\ < 1.0 12]. Surrounding the COT are 
the electromagnetic (EM) and hadronic (HAD) calorime- 
ters, segmented in projective n — cj> towers pointing to 
the nominal collision point z = 0. The central calorime- 
ters measure the energies of particles within the range 
\r\\ < 1.1, while the plug calorimeters extend the range to 
1.1 < I77I < 3.6. Two triggers were used to select the data 
for this analysis. The main trigger requires two high-£ , x 
EM clusters in the calorimeter while a backup trigger ac- 
cepts events with a single electron candidate with very 
high Et and looser electron-selection requirements. 

Events are selected with two high-£r electron candi- 
dates ^| , of which at least one is required to have 
been measured in the central calorimeter. A matching 
COT track is required for all central candidates. Events 
with same-sign central electron pairs are rejected, and an 
isolation condition for the energy found within a cone of 
angular radius R = \J (A7/) 2 + (Acf>) 2 — 0.4 around the 
electron is imposed for electron candidates. The angular 
distribution is measured using cos 9*, where 9* is the an- 
gle between the electron and the incoming quark in the 
Collins-Soper frame (lq. The search is performed using 
events with M ee > 200 GeV/c 2 . 

Z' production is expected to interfere with the SM 
Drell-Yan Z/'y* process, altering the cos#* distribution. 
For this reason, the Z/-f* process is not labeled as a 
"background". Instead, the Z'/Z/j* is referred to as 
the Z 1 signal, while the SM Z/j* is referred to as the 
SM Drell-Yan production. The term background will 
be used to designate all other SM processes (exclud- 
ing Z/j*) expected to contribute to the dielectron final 
state sample. Of these background sources, the most 
important are dijet events in which jets are misiden- 



Source 


Z/-y*^ e+e" 


Dijet 


Diboson 


Total SM 


Observed 


Events 


80.0 ±8.0 


28±i| 


6.8 ± 1.4 
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TABLE I: Summary of the expected and observed numbers of 
dielectron events with M ee > 200 GeV/c 2 in 0.45 fb" 1 . This 
luminosity value was used to normalize the contributions from 
all processes with the exception of the dijet background. 



tified as electrons, and diboson events (see Table QJ. 
The dijet background is estimated using the probabil- 
ity for a jet to be misidentified as an electron ("fake 
rate") which is measured in inclusive jet triggered data. 
The fake rate is then applied to each jet in a sample 
of events containing one electron candidate and one or 
more jets. All non-dijet backgrounds are estimated using 
pythia Monte Carlo simulation, normalized to the 
product of the theoretical cross sections 0, 0] and the 
integrated luminosity. Other background processes such 
as Z/^j* — » t + t~ — » e + e~v T v e u T v e or pair production 
tt — > e + e~ v e i> e bb have negligible contributions in the 
high-mass region considered for this analysis. 

A leading-order calculation is used as the starting point 
to construct the signal and SM Drell-Yan Monte Carlo 
distributions • A next-to-next-to- leading order mass- 
dependent K— factor is then included, followed by a 
(M ee , cos 9*) parameterization of the CDF detector re- 
sponse 20] to dielectron events. This parameterization 
is extracted by running the full detector simulation on 
a large sample of Z/"f* — ► e + e~ events generated with 
pythia such that the distribution in M ee is roughly uni- 
form between 0.03 and 1.05 TeV/c 2 . 

For a particular Z' model (denoted by HI), we iso- 
late the Z' signal by using two variables: the invariant 
mass of the electron-positron pair M ee in 10 GeV/c 2 bins, 
and cos#* in 0.25 bins. The bidimensional distribution 
(M ee , cos 9*) of the CDF data is used to test two mutually 
exclusive hypotheses: 1) the null hypothesis (HO), where 
data are described by SM Drell-Yan and background dis- 
tributions, and 2) the test hypothesis (HI), where data 
points are described by the Z' signal and background 
distributions. A test statistic Q is defined as: 



> p( ^ 1^ \ 

Q(d) = -2- In — W = const 

P(d\H0) 



N bi . 



, KfHl 



where Nun denotes the total number of bins in the 
(M ee ,cos9*) plane, d = (d\, d%, ■ ■ ■ , dN bin ) is the ob- 
served data distribution, while Aff 1 and Mf° are the 
expected numbers of events in bin i in the HI or HO 
hypotheses, respectively. 

Several sources of systematic uncertainty affect our 
measurement. First, a relative uncertainty of 10% on the 
total event rate is incurred primarily due to uncertain- 
ties in the luminosity measurement (6%), the dielectron 
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detector acceptance and electron identification efficiency 
(5%), and the LO calculation (4%). The second dominat- 
ing effect is the electron energy scale and resolution un- 
certainty, which modifies the shape of the M ee and cos 9* 
distributions. The third source is the uncertainty in the 
background (particularly dijet) estimations. The dijet 
prediction uncertainty is extracted from the differences in 
the fake rate measured in kinematically different jet sam- 
ples. Finally, the uncertainty related to the choice of the 
parton distribution functions set (CTEQ6M J2l|) is eval- 
uated using the Hessian method advocated in Ref. [22| . 
and found to have a negligible effect on our results. 

For each of the Z' models (HI test hypotheses) men- 
tioned in the first paragraph, a large number of simu- 
lated experiments d are drawn either from the HI or 
the HO hypotheses, and the corresponding Q(d) values 
are stored in two separate histograms f H1 (Q) and f H0 (Q) , 
respectively. The systematic uncertainties are accounted 
for as described in Ref. j^j. The f H1 (Q) and f H0 (Q) 
distributions are integrated in the region Q > Qq, where 
Qo is the value measured using the d distribution of the 
CDF data. If the ratio of the two integrals is less than 
5%, then the HI model is excluded at 95% confidence 
level (C.L.) H|. 

The CDF data is found to be consistent^ith the null 
(no Z') hypothesis, with a probability P( d \H0) greater 
than the values obtained in 87% of HO simulated exper- 
iments. Figs, n an d El show good agreement between 
data and Monte Carlo SM distributions for the M ee and 
cos#* distributions. For illustration, Fig. 0also presents 
the forward-backward asymmetry A r pg 14] defined as 
(N + - N-)/(N+ +N-), where N + and iV_ are the num- 
bers of forward (cos#*> 0) and backward (cos6>*< 0) 
events in the given M ee range. The superscript "raw" 
is used here to emphasize that no detector acceptance, 
background subtraction, or efficiency corrections are ap- 
plied. The Apb plot is a common way of representing 
the mass dependence of the angular distribution. 

The sequential Z' SM boson, which has the same cou- 
plings to fermions as the SM Z boson, is excluded by our 
data up to a mass of 850 GeV/c 2 , at 95% C.L. It is noted 
here that using the dielectron invariant mass alone would 
require roughly 25% more data for the same Z' SM exclu- 
sion. In general, the improvement provided by including 
the cos 9* spectrum depends strongly on the particular Z' 
under investigation and the integrated luminosity being 
analyzed. Other Z' theories include grand unification Eq 
models 0, 0, , where the E 6 gauge group breaks down 
as £ 6 ->SO(10)xU(l)v, ->SU(5)xU(l) x xU(l) v ,, and the 
SM gauge structure results from breaking down the 
SU(5) group. Therefore, an extra Z' will be a combi- 
nation of the two U(l)'s: Zg Ee = Z^smBse + Z x cos9e b - 
TableHUlists the 95% C.L. lower limits on the M z > for the 
Z x , Zip, Z n , Zj, Zn, and the secluded Z sec Eq models. 

Another class of theories addressing electroweak sym- 
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FIG. 1: M ee distribution of the data (points) compared to the 
prediction for SM Drell-Yan and backgrounds. The individ- 
ual contributions are stacked as follows: other backgrounds 
(dark grey), dijet background (light grey), and SM Drell-Yan 
(open). The inset shows the M ee distribution of high-mass 
data events using a bin size of 10 GeV/c 2 . There are no 
events in the sample with M ee > 500 GeV/c 2 . 



metry breaking and the hierarchy problem are the lit- 
tle Higgs theories 0, where Z' H bosons are predicted in 
order to stabilize the Higgs mass against quadratically 
divergent one-loop corrections. In the minimal model 
of this type (the littlest Higgs), the Z' H couples to left- 
handed fermions only, and these couplings are parameter- 
ized as functions of the mixing angle cotangent cot 9h • 
Our results for cot 9h = 0.3, 0.5, 0.7, and 1.0 are shown 
in Table ITT1 and improve the results reported in 0. 

A recent study reported in defines a more general 
set of Z' models. This study uses simple constraints such 
as generation-independent fermion charges and gauge 
anomaly cancellations to reduce the number of param- 
eters (17) required to define an arbitrary Z' model. Four 
families of models have been considered, each of them 
specified by three parameters: the mass Mz>, the gauge 
coupling g z , and the ratio x of certain U(l) charges. This 
three-dimensional parameter space is sampled to obtain 
the exclusion contours shown in Fig. [3J For small \x\ and 
9z < 0.10, the exclusion limits are more stringent than 
the ones derived in Ref. 0] based on the LEP II results 

Finally, Z 1 mass constraints can be derived from 
searches for contact interactions, if the collider energy 
is far below the Z' pole 0, IH E3- An effective La- 
grangian for the qqee contact interaction can be writ- 
ten as: J2 q T ll , j =L,R 4nr iei^etqjl tJ ,qj/^ J , where A is 
the scale of the interaction, and 77 = ± 1 determines 
the interference structure with the Z/j* amplitudes |26| . 
A generation-universal interaction is assumed and lower 
limits are measured for A in six helicity structure scenar- 
ios: LL, LR, RL, RR, VV and AA (Table III) H3- 

In conclusion, we have searched for Z' decays to e + e~ 
pairs in 0.45 fb _1 of data accumulated with the CDF II 
detector. To strengthen this search, the reconstructed 
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Z' Model 


ZsM 


z x 


z$ 


Z-q 


Zi 




Z EC c 


70.3 


7-0.5 
Z H 






Exp. limit (GeV/c 2 ) 


860 


735 


725 


745 


650 


710 


675 


625 


765 


835 


910 


Obs. limit (GeV/c 2 ) 


850 


740 


725 


745 


650 


710 


680 


625 


760 


830 


900 



TABLE II: Z' exclusion summary: expected and observed 95% C.L. lower limits on M Z ' for the sequential, the canonical E(,, 
and the littlest Higgs Z' models. 
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FIG. 2: Top: Distributions of cos#* for the high mass re- 
gion M ee > 200 GeV/c 2 . The points are the data, the open 
histograms are the predictions from Drell-Yan Monte Carlo 
simulation, and the shaded histograms are the background 
predictions. The individual contributions are stacked. Bot- 
tom: distributions of the forward-backward asymmetry A r pg 
for the data (points) and predicted SM processes (histogram). 
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FIG. 3: Exclusion contours for the B—xL, 10+a;5, d—xu, and 
q + xu Z' families. The dotted lines represent the exclusion 
boundaries derived in Ref. from the LEP II results [24|. 
The region below each curve is excluded by our data at 95% 
C.L. Only models with M z , >200 GeV/c 2 are tested, which 
explains the gap at small for some models. 



Interaction 


LL 


LR 


RL 


RR 


VV 


AA 


A+ limit (TeV/c 2 ) 


3.7 


4.7 


4.5 


3.9 


5.6 


7.8 


A~ limit (TeV/c 2 ) 


5.9 


5.5 


5.8 


5.6 


8.7 


7.8 



TABLE III: 95% C.L. lower limits for the contact interaction 
mass scales. 



dielectron invariant mass M ee spectrum and the angular 
distribution of the electron pair cos 9* are analyzed simul- 
taneously. This is the first study of this kind at the Teva- 
tron, and it opens up a new avenue for exploring the Z' 
production in the inverse femtobarn luminosity regime. 
Many of the Z' models encountered in the literature are 
surveyed, no significant evidence for signal is found, and 
95% C.L. limits are set on these models. Constraints are 
also placed on contact interaction mass scales far above 
the Tevatron energy. Finally, exclusion contours for the 



generic Z' model-lines advocated in Ref. 2] are mapped 
out. In comparison to the LEP contact interaction Z' 
search results given in , our results exhibit higher sen- 
sitivity in the small |x| and small g z regions. 
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